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Abstract 

We present a calculation of the renormalization coefficients of the quark bilinear operators and 
the K — K mixing parameter ■ The coefficients relating the bare lattice operators to those in 
the RI/MOM scheme are computed non-perturbativcly and then matched perturbatively to the 
MS scheme. The coefficients are calculated on the RBC/UKQCD 2+1 flavor dynamical lattice 
configurations. Specifically we use a 16^ x 32 lattice volume, the Iwasaki gauge action at /? = 2.13 
and domain wall fermions with = 16. 
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I. INTRODUCTION 



The RBC and UKQCD collaborations have recently performed the first simulations with 
2+1 flavor domain wall fermions [l, [2I, [sj. Much interesting phenomenology requires the 
conversion of bare lattice quantities to a less arbitrary and more perturbatively amenable 
continuum scheme. In particular, this is true for the determination of weak matrix elements 
such as Bk and for the Standard Model parameters such as quark masses. Of course, 
physical quantities are independent of the choice of renormalization procedure, nevertheless 
theoretical predictions are often given in terms of the parameters of the theory (a^ and 
quark masses) which require renormalization. In addition, for many processes (e.g. K — 
K mixing) the amplitudes are factorized into products of perturbative Wilson coefficient 
functions and operator matrix elements which contain the long-distance effects. The Wilson 
coefficients and operator matrix elements need to be combined with both evaluated in the 
same renormalization scheme. The purpose of this paper is to determine the factors by 
which matrix elements computed in our numerical simulations should be multiplied in order 
to obtain those in the MS scheme which is conventionally used for the evaluation of the 
coefficient functions. 

In principle, for a sufficiently small lattice spacing a and a sufficiently large renormaliza- 
tion scale yU, it is possible to perform the renormalization of the bare lattice operators using 
perturbation theory. However, in practice the coefficients of lattice perturbation theory are 
frequently large leading to a poor convergence of the series and even with attempts such as 
tadpole improvement to resum some of the large contributions, it appears that the typical 
n-loop correction is numerically of 0(q;"), in contrast to continuum perturbation theory 
where the corresponding contributions are of 0((as/47r)"). A related difficulty is the choice 
of the best expansion parameter (a^), for example between some tadpole improved lattice 
coupling or the MS coupling. In practice, at one-loop order, different reasonable choices can 
lead to significantly different results. For the quark bilinear operators and Bk considered in 
this paper, we present the perturbative results and illustrate these points in Section HTl 

The main purpose of this paper is to avoid the uncertainties present when using lattice 
perturbation theory by implementing the Rome- Southampton RI/MOM non-perturbative 
renormalization technique j4|. The key idea of this technique is to define a sufficiently 
simple renormalization condition such that it can be easily imposed on correlation functions 
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in any lattice formulation of QCD, or indeed in any regularization - that is, the condition 
is regularization invariant (RI). We therefore introduce counter-terms for any regularization 
such that a Landau gauge renormalized n-point correlation function with standard MOM 
kinematics at some scale /i^ has its tree level value. This condition is simple to impose 
whenever the renormalized correlation function is known in any regularization. It applies 
equally well to both perturbative expansions to any order and to non-perturbative schemes 
such as the lattice, and thus RI/MOM is a very useful interface for changing schemes. In 
particular, only continuum perturbation theory and the lattice regularization are required 
to obtain physical results from a lattice calculation. 

Our choice of lattice action is important for the efficacy of the RI/MOM technique. With 
domain- wall fermions, 0{a) errors and chiral symmetry violation can be made arbitrarily 
small at fixed lattice spacing by increasing the size of the fifth dimension. This allows us 
to avoid the mixing between operators which transform under different representations of 
the chiral-symmetry group; this is a very significant simplification compared to some other 
formulations of lattice QCD. The action and operators are also automatically 0{a) improved. 

Another important property of DWF is the existence of (non-local) conserved vector and 
axial currents. This will be discussed in detail below. 

In this paper we study the renormalization of the quark bilinear operators ipTip, where F 
is one of the 16 Dirac matrices, and of the AS* = 2 four-quark operator Oil. Table [J contains 
a summary of our results, relating bare operators in the lattice theory with Domain Wall 
Fermions and the Iwasaki gauge action at f3 = 2.13 (a^^ = 1.729(28) GeV, see Section UTTl 
for further details) to those in two continuum renormalization schemes. Columns three 
through five give the three independent Z factors which, when multiplying the appropriate 
bilinear lattice operator, convert that operator into one normalized according to either the 
RI/MOM or MS(NDR) schemes. The final column contains the combination of factors 
needed to convert a lattice result for the parameter Bk into the corresponding RI/MOM or 
MS(NDR) value. 

The plan of the remainder of the paper is as follows. In the following section (Section ln|l 
we start by reviewing the perturbative evaluation of the renormalization constants; the 
results can later be compared with those obtained using the non-perturbative procedures. 
In Section llllt we begin the description of the non-perturbative computations with a brief 
introduction to the details of our simulation and to the computation of the quark propagators 
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which are the basic building blocks for all our subsequent calculations. In Section [IV] we give 
a short introduction to the regularization independent (RI/MOM) scheme. In this section 
we also discuss the renormalization of flavor non-singlet bilinear operators, including the 
check of the Ward-Takahashi identities. The discussion of the renormalization of the four- 
quark operators and the results for the renormalization constant for Bx are presented in 
Section |Vl Section I VI I contains a brief summary and our conclusions. 

II. PERTURBATION THEORY 

Before proceeding to describe our non-perturbative evaluation of the renormalization con- 
stants we briefly review the corresponding (mean field improved) perturbative calculations. 
Specifically, we present perturbative estimates for the renormalization constants of the quark 
bilinears and Bk- These can then be compared to those obtained non-perturbatively below. 
The ingredients for the perturbative calculations and a detailed description of the procedure 

n n 

can be found in refs. p, 16| . 

Writing the domain wall height as M = 1 — cuq, the bare value of Uq in our simulation is 
uq = —0.8. The mean field improved value of cuq is then given by 

to^'^ = uJo + 4:{l-u) - -0.303, (1) 

where the link variable is defined hj u = V^^^ and V = 0.588130692 is the value of the 
plaquette in the chiral limit. 

We define the renormalization constant, Zo^, which relates the bare lattice operator, 
Qhatt^^-i^^ to the corresponding renormahzed one in the MS scheme at a renormalization 
scale of fi = a'^ by: 

Of^{a-') = Z,Of^'\a-'). (2) 

Here i = S,P, V, A, T for the scalar and pseudoscalar densities, vector and axial-vector 
currents and tensor bilinear and i = Bk for the AS = 2 operator which enters into the 
mixing amplitude (or more precisely for the ratio of the AS" = 2 operator and the 
square of the local axial current, which is the relevant combination for the determination of 
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Bk)- The one-loop, mean field improved estimates for the Zi are: 

Zb, = 1-^1.470. (6) 

where Cp is the second Casimir invariant Cp = {N"^ — 1)/2N for the gauge group SU{N). 



Here ^/Z^ is the quantum correction to the normalization factor a/1 — cUq physical 
quark fields (the factors depending on this normalization cancel in the evaluation of Zbh)- 
At one-loop order in perturbation theory 

^. = 1 + ^5.251. (7) 

In obtaining the coefficients in Eqs. (jS]) - ([7j) we have interpolated linearly between the entries 
for M = 1.30 and M = 1.40 in tables III and IV of ref. ^ to the mean-field value of 
M = 1.303. Since the mean-field value of M is so close to the quoted values at M = 1.30, 
we prefer this procedure to using the general interpolation formula quoted in j^. The 
difference between the two procedure is negligible compared to the remaining systematic 
uncertainties. 

In order to estimate the numerical values of the renormalization constants we have to 
make a choice for the expansion parameter, i.e. the coupling constant as- Here we consider 



two of the possible choices, the mean-field value as defined in eq. (62) of ref. and the MS 
coupling, both defined at yU = a~^. The mean field improved coupling constant is given by 

1 V 

9 =-2+dg + Cp + Nfdf , (8) 

where go is the bare lattice coupling constant {qq = Q/f3), and the remaining parameters 
are defined in ref. Gj and take the numerical values dg = 0.1053, Cp = 0.1401 and for 
cjMF ^ _o.303, df = -0.00148. We therefore obtain 

aMF(1.729 GeV) = 0.1769. (9) 

Such a value of the coupling is significantly lower than that in the MS scheme at the same 
scale, for which we take, a^(1.729 GeV) = 0.3138. 



The difference in the two values of the couphng constant leads to a significant uncertainty 
in the estimates of the renormalization constants at this order, as can be seen in Table HTl 
The need to eliminate this large uncertainty is the principle motivation for the use of non- 
perturbative renormalization. The entries in Table [TTl are the factors by which the matrix 
elements of the bare lattice operators should be multiplied in order to obtain those in the 
MS(NDR) scheme at the renormalization scale /i = 1.729 GeV. 

Finally we perform the renormalization group running from fi = 1.729 GeV to obtain 
the normalization constants at other scales, and in particular at the conventional reference 
scale of /i = 2 GeV (see Table [In|) . In each case we use the highest order available for the 
anomalous dimension; two loops for Bx, three loops for the tensor operator and four loop 
for the scalar /pseudoscalar densities. This is the same procedure which we use for the non- 
perturbatively renormalized normalization constants below and the details and references to 
the anomalous dimensions are presented in sections lIV Ft HVGl and IVD] below. The numbers 
in Table UTTl are the factors by which the matrix elements of the bare lattice operators should 
be multiplied in order to obtain those in the MS(NDR) scheme at /i = 2 GeV. The entries 
in the first column indicate which coupling was used in matching between the bare lattice 
operators and the MS(NDR) scheme at = 1.729 GeV, i.e. before the running to other 
scales. 



III. SIMULATION DETAILS 



The calculations described below were performed on the 2+1 flavor dynamical lattice 
configurations generated by the RBC and UKQCD collaborations [2]. The lattices were 
generated with the Iwasaki gauge action at j3 = 2.13 and the domain-wall fermion action 
with Lg = 16. The size of the lattices used in this work is 16^ x 32. The lattice spacing is 
a^^ = 1.729(28)GeV and the residual mass mres = 0.00315(2) in lattice units We have 
3 independent ensembles with light sea quark mass 0.01, 0.02 and 0.03 respectively. The 
strange sea quark mass is fixed at 0.04. For each ensemble, we have used 75 configurations, 
starting from trajectory number 1000 and with trajectory separation 40. 

Following the Rome-Southampton RI/MOM non-perturbative renormalization 
procedure!^, [Sj, the lattices are first fixed in Landau gauge. Then, on each gauge- 
fixed configuration, we measure the point-point quark propagators S{x,Xo) with periodic 
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boundary conditions in space and time, where xq is the source position and x is the sink. 
We have chosen four different sources to generate the propagators, 

xo e {(0, 0, 0, 0) , (4, 4, 4, 8) , (7, 7, 7, 15) , (12, 12, 12, 24)} . (10) 

Next, a discrete Fourier transform is performed on the propagators, 

S (p, Xo) = ^ 5 {x, Xo) exp [-ip ■ {x - xq)] , (11) 

X 

where 

P„ = -j-n^, (12) 

is a four-vector of integers and 

L^ = Ly = L, = lQ Lt = 32. (13) 

For the we take values in the ranges 

G {-2,-1,0,1,2} and G {-4, -3, -2, -1, 0, 1, 2, 3, 4} (14) 

and require that the squared amphtude of the lattice momenta is in the range < < 2.5. 
In this paper, for simplicity of notation we frequently use lattice units for dimensionful 
quantities such as p and m. When we particularly wish to emphasize the nature of the 
discretization errors we explicitly reinstate the lattice spacing, writing for example, {apY or 
{amy. 



IV. RENORMALIZATION OF QUARK BILINEARS 

We now discuss how Green functions computed on the lattice can be used to obtain the 
non-perturbative renormalization constants relating bilinear operators defined on the lattice 
to those normalized first according to the Rl/MOM and then the MS scheme. In the first two 
subsections below. Sections IIV Al and IIV Bl we briefly introduce the definitions and notation 
that we use in the rest of this paper. Some of these are extracted from the earlier RBC 
paper on quenched lattices [Sj] and are included here for completeness. 

The definition of renormalization factors Zg and Zm for the quark wave function and 
mass in the RI/MOM scheme are given in Section IIV A[ The basic amputated quark- 
bilinear vertex functions are defined in Section HV Bl and the conditions defining the RI/MOM 
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scheme are written down. Since our calculations are necessarily performed at finite momenta, 
the effects of chiral symmetry breaking coming from both the non-zero quark masses and 
spontaneous chiral symmetry breaking are visible. We discuss these effects in detail in 
Section IIVCI for the important case of the vector and axial vector vertex functions. 

Next, as a consistency check for our methods, we discuss the accuracy with which our 
off-shell vertex amplitudes satisfy the axial and vector Ward-Takahashi identities in Sections 
IIV DI and llV El respectively. In the later section, the determination of Zs/Zq is also discussed. 
In Section IIVGI we compare the observed scale dependence of with that predicted by 
perturbation theory and interpret the differences as coming from (a/i)^ errors. These are 
removed to determine first Zm^^^'^^^ and then Z^^ . A similar determination of z^^^^^'^ 
and then Z^^^ is presented in Section IIVGI Finally, in Section IIVH| results for the tensor 
vertex renormalization factor Zt are obtained. 



A. Quark mass and wavefunction renormalization 

First, we define the renormalization coefficients for the quark field and the quark mass 
as the ratio between the renormalized quantities and their bare counterparts, 

1 

Qrenix) = Z;^qo{x) (15) 
rriren = Z^TJlQ. (16) 

where q^^n and go are the renormalized and the bare quark wavefunction, and mren and uiq 
are the renormalized and the bare quark mass. With domain-wall fermions, 

niQ = nif + m^es (17) 

where mj is the input quark mass and mj-es is the residual mass. The renormalized propagator 
(in momentum space) is 

^-ren (p, m^en) = ZqSo {p, ^o) |^^^^_/^^ (18) 

where p is the momentum of the quark propagator. 

Since domain wall quarks enter the calculations described here in three different ways we 
must be careful to clearly distinguish their three distinct masses. As described above, our 
calculations are performed with 2+1 flavors of dynamical quarks. We will use the variable mi 
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to label the input mass mj for the light dynamical quarks and for that of the dynamical 
strange quark. Since we often evaluate products of propagators which depend on a third 
quark mass, that mass is labeled m^g^. In some cases the limit rrivai — ^ may be an adequate 
definition of the chiral limit. However, in order to deal with simple results from which a 
weak quark mass dependence has been removed we will often consider the "unitary" case 
f^vai = fni and take the limit mvai = m/ — 0. Of course, in this case remains non- 
zero but since its value is never changed this causes no immediate confusion. Underlying 
the validity of the Rome-Southampton renormalization scheme is the use of infrared-regular 
renormalization kinematics. Therefore, as our renormalization scale /i becomes larger and 
future calculations more accurate, even this weak quark mass dependence will completely 
disappear. 

As discussed in detail in the renormalization technique requires the existence of a 
window of momenta such that 



In practice, however, violation of these restrictions, especially at the higher boundary, has 
to be considered. Due to the spontaneous breaking of the chiral symmetry as illustrated by 
the non-trivial difference between Z^/Za and Z^jZy at low momenta (see section IV C) we 
have to rely on the calculation in the relatively high momentum region, where [a'pf > 1. 
Fortunately, the effects from breaking the restriction imposed by the finite lattice spacing 
a are small and predictable. They introduce an error of O ((ap)^) to the renormalization 
coefficients which can be removed by quadratic fitting to the momentum dependence. A 
more detailed investigation of this issue is presented in j^. 

It is possible in principle to relax the constraint \p\ ^ Aqcd in eq.(19) by performing 
step scaling, i.e. by matching the renormalization conditions successively to finer (and also 
smaller in physical units) lattices. This is beyond the scope of this paper. 

The regularization independent (RI/MOM) scheme is defined such that by adjusting the 
renormalization coefficients Zg and Zm at the renormalization scale /i, and restricting p in 
a suitable window, we have: 



Aqcd < \p\ < a 



-1 



(19) 




(20) 



lim 



1 



Tr {S^^l (P))p2=^2 = 1- 



(21) 



12m, 



ren 
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Imposing these conditions on the lattice and taking into account the dynamical breaking of 
chiral symmetry at low energies, the additive mass renormalization mres and O (a^) lattice 
artifacts, we have the following asymptotic behavior 8|] for relatively large p"^: 

-^Tr {S-l (p)) = p^C^Z, + Zrr^Z, + am,,,} + O {{apf) . (22) 

[ap] 

The left-hand-side of Eg. (1221) at each of the unitary points (m; = mvai) is calculated 
as the inverse of the average over all propagators, where the average is performed over all 
sources and configurations: 

^ Si {p, xo) 



N 

i=l 



n 



source 



(23) 

^0 J ) 

where Usource = 4 and i G {1, 2, ■ ■ • , N} labels each configuration. (For brevity, in this equa- 
tion and in the following we will suppress the subscript "latf'on all the lattice propagators. 
In the following text, unless otherwise specified, all the propagators S {p) without a subscript 
denote the lattice propagators.) Because of possible correlations between propagators with 
different sources calculated on the same configuration, each group of 4 propagators from the 
same configuration is considered as one jackknife bin in the single-elimination procedure. 

In Figure [T] we plot the results for j|^Tr[S'^^(p)] as a function of the momentum and 
tabulate the corresponding numerical values in Table IIV[ In addition to our results for 
non-zero quark mass, we also plot and tabulate the results extrapolated to the chiral limit 
where mi = m^ai = — ?^res for each momentum. Two quantities of interest can be deduced 
from the mass dependence shown in Figure [H First, the chiral limit gives a measure of 
the spontaneous and explicit chiral symmetry breaking and is given in the left-most column 
of Table HVl Second determining the slope with respect to mvai provides one method to 
calculate Z^Zg. This is used later in Section HVEI to test the vector Ward Identity which 
relates this to a second method of computing ZmZ^. 



B. Renormalization of flavor non-singlet fermion bilinears 

We now consider the renormalization of quark bilinear operators of the form uTd, where 
r is one of the 16 Dirac matrices. The corresponding renormalization constant Zr is the 
factor relating the renormalized and bare bilinear operators: 

[urd\,,M = Zr{f^o.)[u^<^o , (24) 
11 



where fi is the renormahzation scale and we will treat only local operators where the lattice 
fields u and d in the bilinear operator [uFdjo are evaluated at the same space-time point. 

Following the Rome-Southampton prescription^ for renormahzation in the RI/MOM 
scheme, we define the bare Green functions between off-shell quark lines, and evaluate their 
momentum-space counterparts Grfi (p) on the lattice, averaged over all sources and gauge 
configurations, 

1 ^ f 1 

. 1 I source ^ 
1=1 \ xo 

We then amputate this Green function using the averaged propagators, 

nr,o (P) = (p) Gr,o (p) (75 [S'' (p)] ^ 75) , (26) 
where 5*"^ (p) is calculated according to Eq. fl23|. The bare vertex amplitudes are obtained 



(25) 



from the amputated Green functions as follows 



A5 (p) 
Ap (P) 

Ay ip) 
Aa (P) 
At (p) 



Y^Tr [Hi (p) 1] 



1 

12 
1 

48 
1 

48 



Tr[n,3(p)75] 
Tr 



Tr 



72 



Yl ^7. (P) 7a< 

(P) 757/. 



- Mis' 



(27) 
(28) 

(29) 
(30) 
(31) 



The values of all the five bare vertex amplitudes at the unitary mass points mi = mvai are 
presented in Table |V] through Table IVIIi Finally, by requiring that the renormalized vertex 
amplitudes satisfy 

1, te{S,P,V,A,T}, (32) 



rpn -i- i-7 



we can calculate the relevant renormahzation constants. 

Equations (12^ through Eq. ( l32l) describe the schematic procedure used to calculate the 
renormahzation coefficients of quark bilinears. In practice however, with finite quark masses 
and a limited range of momenta, we have to consider lattice artefacts and other systematic 
uncertainties. We explain the details in the following sections. 
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C. Chiral symmetry breaking and — Zy 

In this section we examine the effects of both the low-energy spontaneous chiral symmetry 
breaking present in QCD and our non-zero quark masses on the large-momentum, off-shell 
propagators which we are using to impose non-perturbative renormalization conditions. A 
good quantity to study in order to understand these effects is the difference of the off-shell 
vector and axial vector vertex functions. 



1. Numerical results for Za — Zy 

In the limit of a small mass and a large momentum, we expect 

Za = Zy, (33) 

or equivalently, 

Aa{p') = Avip') (34) 

for p2 > Al^cB^ VP?. 

However, with finite quark masses and at relatively low momenta Ky and may receive 
different contributions of the form 

^ (35) 



and 



^Tivai (gg) 



p4 



(36) 



Here we are exploiting the SUi^i^^ x SUii{3) chiral symmetry of large Lg domain wall fermions 
which implies that a difference between Ay and A^ requires the mixing of (8, 1) and (1, 8) 
representations and hence involves a product of two quantities which transform as (3, 3) and 
(3, 3). This requires the two powers of m^ai in Eq. fl35l) and the product m^ai {qq) in Eq. fl36l) . 
The extra factors of and come from naive dimensional analysis. 

To determine how much chiral symmetry breaking is present in our calculation, we exam- 
ine the relative difference between A^ and Ay. In Figure [2] we plot the quantity (jV^^Ay')/2 
as a function of momentum. At relative low momenta, 0.5 < {pa)"^ < 1, we observe that 
this quantity is quite large (~ 5%). Furthermore, even when we extrapolate (-Aa+AvV^ 
the chiral limit, where the terms in Eqs. (135|) and (!36l) both vanish, the difference between 
A A and Ay does not vanish. Here to obtain the chiral limit shown in Figure [2] we perform a 
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linear extrapolation mvai + Tnres 0. While a quadratic extrapolation gives a similar result, 
this linear choice is motivated by the analysis presented in Section IIVC3[ 

Since the explicit chiral symmetry breaking effects needed to split A^i and Ay can be 
argued 9| to be O(m^gg), we would not expect this difference to reffect explicit, finite-L<,, 
domain wall chiral symmetry breaking. In fact, similar deviations between Ay and are 
seen on lattice ensembles without fermion loops where explicit domain wall chiral symmetry 
breaking is expected to be smaller. This is shown in Figure [3] where we plot the same 
quantities from a quenched simulation using the DBW2 gauge action. Thus, it appears that 
this difference represents the high energy tail of QCD dynamical chiral symmetry breaking 
rather than the explicit chiral symmetry breaking coming from the finite value of Lg. 

While the effects of spontaneous chiral symmetry breaking will not vanish in the limit 
'^vai + "^res 0, it is Unlikely that the substantial difference found for A^ — Ay in the chiral 
limit can be explained by a dimension-6 condensate such as 

^ (37) 

since it is suppressed by six powers of momentum and appears to be too small for the size of 
the breaking we have observed. We have also fit the quantity different powers 

of p, as is shown in Figure HI and it is clear that the momentum dependence of the chiral 
symmetry breaking term is dominated by p~'^ or p~^, very different from p"^ that naive 
dimensional analysis suggests should appear in the (gg)^ term above. 



2. Effects of exceptional momenta 

In fact, we believe that the origin of the difference between A^ and Ay is different. 
Our choice of kinematics corresponds to so called "exceptional momenta" , i.e. a momentum 
transfer is zero. This invalidates the naive power counting estimates used above and permits 
the low-energy, spontaneous chiral symmetry breaking to split Ay and A^ with only a 
suppression for large p, as we now explain. Begin by considering a general, amputated 
Feynman graph F with F external fermion lines and B external boson lines. Recall that for 
connected graphs the degree of divergence (i of F is defined as d = 4 — 3F/2 — B. If the 
graph F is disconnected then its degree of divergence is the sum of those of its connected 
components. Now imagine that each external line of F carries an incoming momentum A Pi 
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for 1 < i < F + B , where A is an over-all scale factor. The asymptotic behavior for large 
A of the amplitude corresponding to such a graph will be A'^' where d' is the degree of 
divergence of a subgraph F' C F. This subgraph F' must be chosen so that i) there exists a 
routing of the internal momenta within F such that all lines carrying momenta proportional 
to A lie within F' and ii) F^ possesses the least negative degree of divergence d' of all those 



subgraphs satisfying i) [lO|, . Note, that F' may equal the original graph F and may itself 
be disconnected. 

The most familiar situation is the case of non-exceptional momenta, defined as a mo- 
mentum configuration in which no proper partial sum of the external momenta pi vanishes. 
Under these circumstances all subgraphs F' obeying i) must be connected. (Otherwise there 
would be zero momentum transfer between the groups of momenta entering each of the dis- 
connected components.) This implies that the subgraph F' with the least negative degree of 
divergence is one with no additional external lines beyond those already appearing in F which 
in turn implies that this subgraph F' is the entire graph F. For the case of the vertex graph 
of interest, we deduce a constant behavior (up to logarithms) since (i = 4 — 1 — 2-3/2 = 0. 
(Here it is convenient to view this vertex graph as resulting from a normal Feynman graph in 
which an external vector boson is coupled to the vertex so the rules discussed above directly 
apply.) 

This analysis not only gives the leading asymptotic behavior but also insures that extract- 
ing a few extra factors of the mass m or the chiral condensate (qq) will make the degree of 
divergence of that graph more negative and hence make its asymptotic fall-off more rapid, in 
the fashion suggested by naive power counting. For the case of interest, we would like to re- 
strict a subset of the internal fermion lines of our graph F to carry only low momenta so that 
they will reflect the low-energy, spontaneous chiral symmetry breaking of QCD. By defini- 
tion, these low momentum lines cannot enter the subgraph F' discussed above whose degree 
of divergence determines the asymptotic behavior of amplitude being studied. In order to 
split Ay and A^, chiral symmetry breaking transforming as an (8,8) under SU{3) x SU{3) 
is required. This in turn requires that this low energy, excluded subgraph must be joined to 
the remainder of the graph by at least four fermion lines. 

Such a circumstance is illustrated by the general vertex graph F in Figure [5l contained 
in the outer dashed box. Here we have identified a subgraph F2 which carries only low 
momenta and can therefore transform as (8,8) even in the limit of vanishing quark mass. 
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^vai+'^res = 0. Foi the case of non-exceptional momenta, we must apply Weinberg's theorem 
to the subgraph F', enclosed in the inner dashed box, through which, by assumption, all of 
the large momenta entering the vertex and the two external fermion lines must be routed. 
Because of its connections to the subgraph the subgraph T' has six external fermion lines 
and one external boson line (connected to the vertex). The resulting degree of divergence is 
(i' = 4 — 1 — 6-3/2 = —6, justifying the naive behavior in Eq. [371 

However, in our case Ay and are being evaluated with zero momentum entering 
the current vertex and with a vanishing sum of the two incoming fermion momenta — a 
configuration of exceptional momentum. For such a choice of external momenta we can 
divide the subgraph F' identified above into two pieces Fi and F3. Because the momenta 
are exceptional with no large momenta entering the vertex, we can route all of the large 
momenta through F3. Since F3 has only four external fermion lines, its degree of divergence 
is = 4 — 4-3/2 = —2 and the behavior above has been replaced by the much less 
suppressed l/p"^. If we think of the subgraph F2 as a generalized chiral condensate {0\qqqq\0) 
we are seeing the asymptotic behavior 

very consistent with our numerical results. Note the discrepancy in dimensions between 
Eqs. ( 1371) and ( l38l) will be made up by four powers of Aqcd, the momentum scale to which 
the subgraph F2 is restricted. 

A simple class of graph allowing this behavior can be seen in Figure l6l Here the large 
momentum carried by the two external fermion lines can be routed through the gluon propa- 
gator that is shown explicitly so that the upper part of the diagram carries only low momenta. 
The large momentum behavior of the gluon propagator gives the expected 1/p^ behavior. 
The two general fermion propagators shown with the shaded "blobs" carry small momenta 
and, as suggested by Eq. (122|) . can show (3, 3) or (3, 3) chiral symmetry violation even when 

To confirm this analysis, we have also calculated the difference between A^ and Ay with 
non-exceptional momenta. We have chosen 5 different momentum scales, each corresponding 
to a set of momenta which satisfy the condition pi = p^ = (pi — ^2)^ = for five values 
of p^, as listed in Table IVIIII and Table IIXI We then calculated A a and Ay with the two 
external fermions carrying respectively pi and p2- The result is plotted in Figure 13 which 
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shows that the chiral symmetry breaking vanishes almost completely with non-exceptional 
kinematics at medium to large momenta. 

While it would be more satisfactory to perform the calculations presented in this paper 
using non-exceptional momenta, the resulting RI/MOM normalization conditions would not 
correspond to those for which perturbative matching calculations have been carried out. 
Thus, we would not be able to relate the quantities which we calculated to those defined in 
the MS scheme. (Of course, this difficulty will be removed when the necessary perturbative 
calculations have been extended to non-exceptional kinematics.) A second, less significant 
advantage of the exceptional momenta which we use is that the exceptional momentum 
conditions are satisfied by a much larger set of discrete lattice momenta permitting the 
RI/MOM condition to be satisfied for more fine-grained sequence of energy scales. 

We now return to the calculation with exceptional momenta {pi = P2), at the scale which 
we are most interested in, that is ~ 2 GeV or {apY ~ 1.3, where and Ay have a 
difference of about 1%. Since we have no means to determine which of these two quantities 
has less contamination from low energy chiral symmetry breaking we have decided to take 
the average | (A^ + Ay) as the central value for both Zg/ZA and Zg/Zy. The difference 
between Aa or Ay and | (A^ + Ay) then provides an estimate for one systematic error in 
our final results. The value of ^ (A^ + Ay) is plotted in Figure [H 

3. Chiral extrapolation to vanishing quark mass 

As discussed above, our use of exceptional momenta implies a 1/p^ suppression for both 
terms behaving as mvai and ^^^1- '^^^ added dimension of a m^^^j term can be provided by 
a factor of 1/Aqcd without the need to introduce additional inverse powers of p. For our 
largest value of m^aia = 0.03, we might estimate mvai/AgcD ~ 0.2. This suggests that we 
should expect a linear rather than quadratic behavior in rrivai to dominate the small quark 
mass limit. 

The difference A a — Ay discussed above provides a good place to study this effect. This 
difference reflects the chiral symmetry breaking of interest and may make these effects stand 
out with possibly reduced errors because of the statistical correlations between the two 
quantities being subtracted. Figure M compares linear and quadratic flts to the dependence 
on the quark mass evaluated at unitary points with m-vai = m; + mres for p = 2.04 GeV. In 
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Table |X] we present the results of these two fits: 



(A.H-Av)/2 + (39) 

(A. + A.)/2 ==» + ^=7^- 

for Aqcd = 319.5 MeV. As can be seen in the Table the linear fits are favored. The linear 
fits have the smaller and the coefficient ci is significantly closer to an expected value of 1 
than is the coefficient C2- Thus, based on both the theoretical expectation and this empirical 
evidence, we will adopt this linear description in the remainder of this paper and extrapolate 
our exceptional momentum amplitudes to the chiral limit using a linear ansatz. For the case 
at hand. Figure [TO] shows this linear extrapolation for the average | (A^ + Ay) to the chiral 
limit for the momentum p = 2 GeV. Figure [H] shows the results in the chiral limit as a 
function of momentum. The results in the chiral limit are also presented in Table IXIi 

D. Axial Ward-Takahashi identity 

Performing an axial rotation on the propagator leads to a relation between Ap and 
Tr (5*^^), the axial Ward-Takahashi identity|8|]: 

A,(„.1JL™_. (41) 

12 (m^ai + mres) 

In a truly chiral theory or the present DWF calculation in the limit — > oo (when chiral 
symmetry becomes exact and m^es = 0), this identity will be obeyed configuration by config- 
uration. However, for finite Lg and mres 7^ 0, this relation will hold only after a gauge field 
average, {e.g. m^^s is only defined after such an average). Thus, we should check Eq. (14T|) 
on gauge-averaged amplitudes. 

Figure [TT] shows the difference between the l.h.s and r.h.s of Eq. fl4ip . divided by their 
average. The case of mi = 0.01 shows relatively larger breaking of < 8%, while the other 
two masses result in a smaller breaking. Since for m = 0.01, the mres term, with a value of 
0.00308, represents a 30% effect, this suggests that the use of mres in the context of Eq. (14T|) 
may be accurate at the 25% level for this lattice spacing. Note, we expect violations coming 
from a dimension-five, anomalous chromo-magnetic term to be suppressed by a factor of 
(pa)^ relative to those from mres, making this < 8% estimate comfortably smaller than the 
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naive estimate of (mres/O.Ol) ■ (pa)^ ~ 30%. However, the suggested growth in the size of 
these violations with increasing (pa)^ may be visible in Figure [TTl 



E. Vector Ward-Takahashi identity and the chiral limit of A5 

Similar to the case with Ap, from the continuum vector Ward-Takahashi identity, we 

_ n 

have the relation between and Tr {S [8!]: 

A, . 2_0J,flM. (42) 

12 Sm^ai 

We are able to check our data against this identity using the three sources, (0,0,0,0), (4,4,4,8) 
and (12,12,12,24) since it is only for these three sources that multiple valence mass data are 
available for each sea quark mass. 

Figure [12] shows the difference between the two sides of Eq. (H2|) divided by their average. 
For all three sea quark masses the data agrees well with the vector Ward-Takahashi Identity 
(Eq. fH2]) ) for medium to large momenta. Equation (1421) implies the relation 

Zm = ^, (43) 

and will use this equation and a calculation of A5 to determine the mass renormalization 
factor in the following sections. 

To extrapolate A^ to the chiral limit, we will improve upon the discussion in 8(| in two 
regards. First, as explained above, we will exploit the asymptotic properties of Feynman 
amplitudes evaluated at exceptional momenta and assume that the leading mass dependence 
in the chiral limit will be linear in m. This is different from the dependence assumed in 
Ref. jsl where dimensional arguments, appropriate to the non-exceptional case and leading 
to the behavior in Eq. [35] were adopted. 

Second, in contrast to that earlier quenched calculation we can examine the behavior 

of Zs as a function of both the valence and light dynamical quark masses, m^ai and mi 

respectively. In Figure [13] we plot As as a function of both mvai and mi. The three curves 

are each a linear plus double pole fit to the valence quark mass dependence of the form: 

2 

TTl 

As{m^i,h m) = Co{mi) + ci{mi)m^^i + Cdp^^, (44) 

where we have allowed the coefficients cq and ci of the constant and linear terms to vary with 
the dynamical light quark mass. However, we have used a common double-pole term with 
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the mf behavior expected for a theory with two hght flavors. Recall that this double pole 
term arises from topological near-zero modes which for two light flavors will be suppressed 
by two powers of the light quark mass. The data in Figure [T^ shows just this behavior with 
the sharp turn-over at the smallest value of mvai increasing as the light dynamical mass mi 
increases. 

This double-pole can be deduced from Eq. |42l As discussed in Ref. js], the NLO, l/p"^ 
term derived from an operator product expansion of the quark propagator on the right- 
hand side of this equation is proportional to the chiral condensate {qq) Q • Isolated, 
topological, near-zero modes of the sort that arise from a gauge fleld background with 
non-zero Pontryagin index will contribute a term to the chiral condensate which behaves 
as l/mvai- This implies that the derivative in Eq. will yield the double pole, 
hypothesized in Eq. HH Such a near-zero mode will also introduce a factor of mi into the 
fermion determinant of the QCD measure for each light flavor in the theory, hence the 
expected factor of mf in the numerator of Eq. HH In Figure [H] we show the variation of the 
double pole coefficient with the momentum at which the coefficient of the double pole was 
extracted. Also shown in this figure is a fit to the expected l/p^ behavior which describes 
the results very well. 

This understanding of the double pole terms suggests that a good strategy for extracting 
the chiral limit of As first takes the limit of vanishing mi to remove this NLO double pole 
term and then extrapolates to m^ai = 0. In the present case, we perform the simpler linear fit 
using the unitary points to obtain As = Zq/Zs since we do not have the complete partially 
quenched results for each of our four sources. 



F. Mass renormalization and renormalization group running 

To calculate the mass renormalization constant Z^, as defined in Eq. flTBl) . we can either 
directly take the derivative of Tr [S'i~j\ (p)] , by following Eq. fl22|) . 



.„.,.1MJ (45) 

12 am^ai 



or we can use the Ward-Takahashi identity. 



Zm = (46) 
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With the analysis described in the above sections, we find that the method with the smallest 
statistical uncertainty is to use l/Z^ as the value of Zm- To remove the factor Zq from Zq/Zs 
(which is equal to A5 and can be calculated as described in Section nVBp . we use the ratio 
Zq/Z^ calculated in Section [IV CI as well as the value Za = 0.7161(1) obtained in Ref. 
using hadronic matrix elements. We therefore determine Z^ in the RI/MOM scheme by 
computing separately the three factors on the right-hand side of 



^RI/MOM 



|^(P) 



1 



(47) 



Table IXIII contains the values of Z^^^'^^ (p) for a variety of momentum scales. 



After obtaining the lattice value of Z^^^'~'^ at different momenta, we divide it by the 



predicted renormalization group running factor to calculate the scale invariant quantity Z^. 



The four- loop running formula we use is 



3: 



c{as (/u) /tt) 

where /zq is chosen such that (a/io)^ = 2, a value that lies within the fitting range used 
below. For completeness we present in the appendices the detailed procedure for running 
as at four-loops (Appendix \^ and the form of running factors (Appendix [B]) . 

As Figure [T5] shows, the quantity Z^ is remarkably independent of the scale /i. However, 
in spite of the name, for other cases, the scale-invariant Z factors do show noticeable scale 
dependence and an additional correction is warranted. (See, for example. Figure [T71) We 
believe that the primary reason for this lack of scale invariance is the presence of lattice 
artifacts, namely the finite lattice spacing which introduces a small error of O (^{afi)'^Y Such 
an error can be reduced by removing the fi"^ dependence in Z^. To do this we fit this 
momentum dependent Z^ to the form A + B (ajj,)'^ over the momentum range 1.3 < (a/i)^ < 
2.5 and then take the (a/i)^ —>■ limit of that fit to remove the /i^ momentum dependence. We 
interpret the outcome as the true Z^. Note, we are ignoring possible fi dependence arising 
from the absence of higher order terms in the matching factor. Such scale dependence can 
only be removed by even higher order computation of the perturbative matching factor and 
such a correction is expected to be very small. While this procedure represents a negligible 
correction for this case of Z^, it will have a more significant effect in the cases considered 
below. 
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Our ultimate goal is to determine which connects the bare lattice quark mass to its 
continuum counterpart defined according to the MS scheme, at the renormalization scale 
fi = 2 GeV, because the corresponding continuum renormalization is conventionally done 
in this scheme. So we again use Eq. (l48l) to calculate Zm^^'~'^{2 GeV) from the scale- 
independent value of Z^. Then we multiply it with the three-loop matching factor, which 
will also be explained in Appendix [HI to match the Z^^^'~'^ {2 GeV) to the MS scheme. The 
final step is shown in Figure [16] and the results are given in Table Kill The renormalization 
constant at the desired scale is 

Z^ (2 GeV) = 1.656 ± 0.048 (stat) ± 0.150 (sys) . (49) 

The systematic error is determined by adding in quadrature our estimates of three different 
types of systematic error which we will now discuss. 

The first is the effect on Z^ of the difference between determining Zg/Z^ from 
1 (A^ + Ay) or from A^. This contributes an error of ±0.011 to Z^^ (2 GeV). Next we 
must assign an error to our use of three-loop matching factor, given in Eq. flB6l) . Here we 
assign an error equal to the magnitude of the final, order error in this perturbative ex- 
pression, which is ±0.103. While this may be a conservative estimate of the omitted terms 
of order and higher, it also is intended to include the errors introduced by the order 
estimate of the perturbative running determine the intermediate SI step used to remove the 
(a/i)^ errors. 

Finally we address the errors arising from our failure to extrapolate to the limit of van- 
ishing strange quark mass. Recall, we have evaluated the chiral limit in which both the 
valence quark mass which enters our off-shell propagators and the dynamical light quark 
mass are extrapolated to zero. However, all of the gauge ensembles used in this calculation 
were computed with a non-zero strange quark mass = 0.04. Since we are matching 
our Green functions to those computed in perturbation theory in the mass- independent, 
m 0, limit our non-zero value for implies an additional systematic error. Because 
the dynamical quarks enter only through loops, their effect is different from that of the 
valence quarks discussed above. They do not contribute chiral symmetry breaking effects 
in our matrix elements. However, because of low energy chiral symmetry breaking, we do 
expect the dynamical quark masses to appear linearly in a quantity such cLS Zjji ill the limit 
m ^ 0. To estimate the size of this 0{ms) effect, we begin with the size of the observed 
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linear dependence, dZm/dm ^ 5.4 which comes from both the calculated valence and light 
dynamical mass dependence of A5. This is then multiplied by 1/2 because there is only 
one flavor of strange quark and by = 0.04 giving an error in Zm of ±0.108. The total 
systematic error given in Eq. HHl ±0.150, is then the sum of these three errors in quadrature. 



G. Quark wavefunction renormalization and renormalization group running 

In Section llVCt we calculated the ratio of renormalization constants 

1^ = ^ (Aa ± Ay) . (50) 

To calculate Zg, we multiply this quantity with Za = 0.7161(1) obtained in Ref. [7|. Thus, 
we have evaluated the quantity Zg in the RI/MOM scheme, which is shown in Table IXlTTl 
To calculate Zg in the MS scheme, we follow a similar procedure as in the previous section, 
and start by dividing ^^^/^^^^ i^y ^j^g perturbative running factor. As shown in Appendix [Cl 
the functional form of the running factor is quite similar to that of Z^- The energy scale /io 
where Z^^ is fixed to the value is again chosen such that {a^of = 2. 

The calculated values of Z^^ vary slightly with momentum due to the presence of lattice 
artifacts. To remove these, we again fit the dependence to the form A + B {aiif and extrap- 
olate to a = 0. The procedure is shown in Figure [T71 Finally, we take the scale-invariant 
Zg^, run up to different scales in the RI/MOM scheme, and then apply the perturbative 
matching factor (Appendix [C|) to translate it to the MS scheme. The MS values are shown 
in Figure [m and Table IXlTTi Of particular interest, the value at /i = 2 GeV is 

Zf^ (2 GeV) = 0.7726 ± 0.0030 (stat) ± 0.0083 (sys) (51) 

The systematic errors are estimated using the same procedure explained in Section lIVFi 
They are the sum in quadrature of the estimated errors arising from the difference A^ — Ky 
(0.0061), the use of a perturbative matching factor accurate to order (0.0045) and our 
use of a non-zero sea quark mass (0.0035). 



H. Tensor Current Renormalization and Renormalization Group Running 

To calculate the tensor current renormalization constant Zt-, we follow a procedure similar 
to those of the previous two sections. For each dynamical quark mass, we combine the ratios 
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Zq/Zx and Zq/ZA in order to obtain the ratio of Zt to Za in the RI/MOM scheme: 

Zj' 



^RI/MOM 



Za 

Ultimately, we use the independent hadronic matrix element calculation of Z^, which gives 
Za = 0.7161(1), to obtain Z^- Table IXTVl shows the values obtained for the ^^^^^^ the 
chiral limit for a range of lattice momenta. As discussed above we have performed the chiral 
extrapolation using a linear functional form, and Figure [19] shows this linear extrapolation 
at the lattice momentum (ap)^ = 1.388. 

We obtain SI values for Zt in the chiral limit by dividing out the tensor current per- 
turbative running factor, the evaluation of which is described in Appendix [El Again, the 
SI values obtained in this way exhibit a dependence on the lattice momenta, and again we 
fit the momentum-dependent Zj^ to the form A + B[anY and extrapolate to (a/i)^ 
to remove the lattice artifacts, as shown in Figure [201 Finally, we run the scale-invariant 
Zx/Za back to different scales in the RI/MOM scheme and use the perturbative matching 
factor (Appendix [E| to match to the MS scheme. The MS values are shown in Figure [211 
and Table IXIVI At yU = 2GeV, we obtain: 

Z^(2GeV) = 0.7950 ± 0.0034(stat) ± 0.0150(sys). 

The systematic errors are determined in the same fashion as in the previous two sections. 
Specifically the errors arising from the difference A^ — Ay (0.0054), the use of a perturbative 
matching factor accurate to order a (0.014) and our use of a non-zero sea quark mass (0.0003) 
are added in quadrature. 



V. RENORMALIZATION COEFFICIENTS FOR Br 

A. General procedure for computing the mixing coefficients 

By the renormalization of we mean the calculation of the renormalization coefficient 
for the operator 

Ovv+AA = {s^^ (1 - 75) d) (s7^ (1 - 75) d) (52) 

which is the operator responsible for the mixing between and K'^. Since for finite Lg 
our theory does not posses exact chiral symmetry we must consider the possibility that this 
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operator can mix with the four other AS = 2 operators with a different chiral structure: 

Ovv-AA = (sY (1 - 75) d) (s7^ (1 + 75) d) (53) 

Oss+pp ={s{l- 75) d) {s (1 - 75) d) (54) 

Oss-pp = (^ (1 - 75) d) {s (1 + 75) d) (55) 

Ott = (sa^'d) (sa^.d) (56) 

where they are labeled by the chirality structure of the even-parity components. The odd- 
parity components of these operators are not important here since they don't contribute to 
^ mixing. 

For domain-wall fermions, the mixing of Ow+aa with these four operators with wrong 
chirality should be strongly suppressed by 0{'m'^^^). However, chiral perturbation theory 
predicts that the B parameters of the operators with the wrong chirality diverge in the chiral 
limit Q. To address this issue, we will describe a theoretical argument to estimate 
the size of the mixing terms and an actual calculation of these chirality-violating mixing 
coefficients on the 2+1 flavor dynamical lattices. 

Following the Rome-Southampton prescription , we first calculate the 5x5 matrix, 

BADC 

I 

(iaS'f 

where Fj^** is the amputated, four-point Green function. The Green functions are first aver- 
aged over all sources and configurations, and then amputated using the averaged propagator, 
in a procedure similar to the calculation of two-point amputated Green functions nr,o {p) in 
Section [iVBI Pj is a suitable projector, which projects out the component with the expected 
chirality (for example, the projector corresponding to Ow+aa is 7^ ® 7^ + 7'^'75 ® 7/^75)- 
The subscripts z, j G {VV + AA, VV - AA, SS - PP, SS + PP, TT}. 

It is straight-forward to calculate the mixing matrix at tree level which we denote as: 

F., = A [rri . (58) 

The RI/MOM renormalization condition which we adopt is then: 

■^ZijM.k = Fik (59) 



= p, [F-] = (F-)^;;;" (p,)27' (57) 



or 



i-Z = FM-\ (60) 
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B. Theoretical argument for the suppression of mixing coefficients 



As can be seen from the structure of the four operators in Eqs. (153|) . (l55l) and (l56l) . 
if they are to mix with Ow+aa defined in Eq. flS^ then two quark fields much change 
chirahty from left- to right-handed. For domain wall fermions such a mixing can arise from 
the explicit breaking of chiral symmetry coming from the finite separation between the left 
and right walls. The asymptotic behavior for large of the resulting mixing coefficients 
can be estimated using the transfer matrix T for propagation in the s-direction introduced 
by Furman and Shamir il7i|. The large-L^ limit is then controlled by matrix elements of the 
operator T^" and is dominated by those four- dimensional fermion modes corresponding to 
eigenvalues of the transfer matrix which lie near unity. 

As described in detail in Ref. [l^ and in the original references cited therein, these fermion 
modes are believed to fall into two classes: modes localized in space-time with corresponding 
T eigenvalues falling arbitrarily close to unity and de-localized modes characterized by a 
mobility edge Ac > and with eigenvalues of T lying below e"^*" 19|, l20|, l2l|, |22|. Since 



two quarks must change chirality to produce the required operator mixing, for the case of 
de-localized modes such mixing will be suppressed by the two factors of e''^"^" needed for 
the propagation between the left and right walls of these two fermions, consistent with our 
estimate above that such effects should be of order m^^g. 

However, the effects of the localized modes are more subtle. We must address the pos- 
sibility raised by Golterman and Shamir 23|] that the contribution of such a mode to mres 
is suppressed because such modes are relatively rare and the necessary coincidence with the 
location of the operators being mixed is unlikely. However, if present, such a mode can mix 
right- and left-handed fermions with little further suppression since the corresponding T 
eigenvalue may be very close to unity. This raises the possibility that a single such mode, 
suppressed by a single factor of rrires might be occupied by the two different quark flavors 
to provide the double chirality flip needed to mix the operators. Fortunately, as argued in 
Refs. and [ist, this is not possible because the mixing in question requires both a quark 
and an anti-quark or two quarks of the same flavor to propagate across the fifth dimension. 
This requires two distinct modes and hence incurs the double suppression which is well 
represented by the m^^g estimate above. Note, m^^g ~ 10^^ which will introduce 0(0.1%) 
errors in current calculations of Bk 2^ and will be too small to be seen in non-perturbative 
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studies presented here. 



C. Lattice calculation of mixing coefficients 

With the procedure described in Section IV At we can directly calculate the mixing co- 
efficients. In particular, we have calculated the off-diagonal terms in the matrix FM~^. 
Figure [22] shows the mixing coefficient FMyy^j^^ vv-aa different unitary masses. As in 
the earlier discussion of the A^ — Ay difference, our use of exceptional external momenta 
permits both a linear and quadratic mass dependence. As was found in Ref. [l^ and sug- 
gested by the mass dependence seen in Figure [23] a linear dependence appears reasonable 
and it is a linear form that we have used in determining the chiral limit shown in Figure [22] 

As can be seen in Figure [22] at the chosen reference scale, fi 2 GeV or (ap)^ ~ 1.4, 
the mixing coefficient is about 0.7% and decreases when the scale is made larger. Similar to 
the discussion in Section flV CI we again propose that this non-zero mixing coefficient in the 
chiral limit has its source in our use of exceptional momenta. Again we can determine the 
asymptotic behavior of the amplitude in question by determining the least negative degree 
of divergence of a subgraph F' through which all of the large external momenta can be 
arranged to flow. Here it is convenient to treat the operator Oll, which is evaluated at zero 
momentum, as an internal vertex of dimension 6 rather than an unusual sort of external 
line. This alters the rules for computing the degree of divergence of a subgraph: now any 
connected subgraph with F external fermion lines and B external boson lines in which this 
new Oil vertex appears, must have degree of divergence d = 6 — 3F/2 — B since On has 
a dimension two higher than the usual renormalizable coupling. (As before, the degree of 
divergence of a disconnected graph is the sum of the degrees of divergence of its connected 
components.) 

As in the case of the vertex amplitude discussed in Section IIVCI the appearance of 
exceptional momenta does not change the asymptotic behavior in the large A limit with 
external momenta \pi for 1 < i < 4. Even for our exceptional case Pi = Ps = —p2 = —P4, 
X'^ scaling with (i = 6 — 4-3/2 = is expected. However, derivatives with respect to the 
quark mass or the occurrence of factors of (qq) will be strongly affected by this choice of 
external momenta. As is shown in Figure [211 "we can identify a disconnected subgraph F' 
through which all the large external momenta can be routed which has (i = (6 — 4- 3/2) + 
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(4 — 4-3/2) = —2. (Note, momentum conservation implies that if all of the large momenta 
can be routed within a disconnected diagram then the choice of external momenta must be 
exceptional.) This d = —2 value implies a l/p"^ behavior with only low momenta flowing 
through the omitted subgraph Fi. Since Fi has four external lines it can translate standard 
QCD vacuum symmetry breaking into the chiral symmetry breaking that is required to 
produce the operator mixing shown in Figure [221 

Again, we confirm this conclusion, by recomputing the coefficient -^^yy+AA.vv'-AA 
non-exceptional momenta, as shown in Figure [251 With that choice of momenta the mixing 
coefficient vanishes completely within our statistical accuracy. 

The other chiral symmetry breaking mixing coefficients, FMyy_^_^^ ss±pp 
FMyyj^j^j^ rprp 8.16 veTj slmllar to the case of FMyy_^_^^ vv-aa J^^^ discussed. These co- 
efficients are plotted in Figure [26] to Figure [291 Since our theoretical argument implies that 
the mixing coefficients are very small, i.e. O {fn^^J and our numerical results are consistent 
with this implication, it is safe to neglect them and calculate the renormalization coefficient 
for Bk: 

^•RI/MOM _ Zvv^AA,VV+AA /g-j^x 
D. Calculation of Zbj^ and renormalization group running 

Using Eq. ( l6Tll . the value of Zg/ZA from Section IIV CI and the value of 
Zg"^ Zvv+AA,vv+AA = F^vv+AA vv+AAi ^^"^ Calculate the lattice values of Zbj^ at dif- 
ferent masses and momenta, as shown in Table [XVl To extrapolate to the chiral limit, we 
again use a linear function, for the same reasons as described in Section IIV CI The linear 
mass fit at the scale ji = 2 GeV is illustrated in Figure [301 and the value of ^'^^^^^ in the 
chiral limit is shown in Figure [3ll and Table IXVII 

Similar to the procedure described in Section IIVFI in order to determine from 
Z^^^^^ , we first divide the Z^j^^°^(/i) by the predicted running factor at one-loop order 
and obtain the quantity Z^^{^). Then we fit a quadratic function A + B (afi)'^ over the 
region 1.3 < (a/x)^ < 2.5 to remove the O [{afi)'^) dependence from Z^^ induced by the 
lattice artifacts. Finally, we restore its perturbative running in the MS scheme to the scale 
fi = 2 GeV. The perturbative running and matching factors are presented in Appendix [Dl 

The procedure of dividing by the running and removing the (a/i)^ dependence is shown 
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in Figure 1321 and the result of restoring the running in the MS scheme is shown in Figure 1331 
Table IXVll hsts at different momentum scales. The final Zbj^ we need in the MS scheme 
and 11 = 2 GeV is 

(2 GeV) = 0.9276 ± 0.0052(stat) ± 0.0220(sys). (62) 

The systematic error is calculated, following the same procedure as has been used for the pre- 
vious quantities, as a sum in quadrature of the amount the result changes when ^ (A^ + Ay) 
is replaced by A^ (0.0131), the size of the highest order perturbative correction being made 
(here of 0{as)) (0.0177) and the effect of our non-zero value for m<j in the calculation of 
Zb^ (0.0007). 



VI. CONCLUSIONS 



We have presented a study of the renormalization coefficients Zq, Z^ Z^ and on the 
16^ x32, 2+1 flavor dynamical domain-wall fermion lattices with Iwasaki gauge action of /? = 
2.13 and a'^ = 1.729(28) GeV generated by the RBC and UKQCD collaborations. These 
coefficients are important components in calculations of a number of important physical 
quantities reported elsewhere 3|, |2J] The procedure closely follows that used in an earlier 



study with quenched lattice configurations [8|, ll5||. In addition to providing the Z-factors 
necessary to support a variety of calculations on these lattice configurations, this paper also 
presents a number of new results which go beyond earlier work. 

First, the troublesome double pole which appears in a quenched calculation of the quan- 
tity As because of topological near zero modes is now highly suppressed by the 2+1 flavor 
determinant. This allows us to use As for an accurate calculation of Zm- Second we have 
identified the 0(5%) large chiral symmetry breaking effects seen in the off-shell Green func- 
tions Ay and A^ as caused by our use of exceptional momenta. We have advanced both a 
theoretical discussion explaining the pattern of symmetry breaking which we have observed 
and a calculation with non-exceptional momenta in which these effects are dramatically 
reduced. 

Third, for Zb^^ we have presented both a theoretical argument and numerical calculations 
showing the mixing coefficients with the operators with the wrong chirality are very small so 
that the calculation of Zbj^^ can be simplified by neglecting these mixing coefficients. Finally 
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we have exploited the earher perturbative work of others and evaluated the factors relating 
the normalization of operators defined in the MS and RI/MOM schemes determining Z^^, 
Z^^, Z^^ and Z^^ from their non-perturbative RI/MOM counterparts to three, three, one 
and two loops respectively. 
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Appendix A: THE QCD f3 FUNCTIONS AND THE RUNNING OF a. 



The four-loop QCD beta functions is calculated in 25|] and the conventions we use are 
the same as in 14l |: 



/J.=i(ll-f.v). 

ft 4 (102 

1 /2857 5033 325 . 
f3^= I . , .2 



: Uf H Uf 

64 V 2 18 ^ 54 ^ 



/?3 
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256 



149753 
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3564C3 



1078361 6508, , 
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/ 50065 6472 
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162 ■ 81 'V ^ 729 



162 
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(Al) 



To calculate the coupling constant (/i) at any scales, we have used the four-loop (NNNLO) 

I I 

running formula for as [25j]: 

da. 



din /i^ 



(3 (as) 

-(3oal - (3^al - P^a^, - ^^a^ + O (a^) 



(A2) 
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where = ols/t^- (We have changed the normahzation of to match the definition of the 
/^-functions coefficients.) For a numerical implementation, we start from the world- average 



value at /i = Mz [26 1 



{Mz) = 0.1176 ± 0.002, (A3) 

where the superscript indicates that it is in the 5-fiavor region, and run as across the nib 
and rric threshold with the matching conditions: 

(nib) = a'f^ (m^) and a^J^^ (m,) = ^ (m,) . (A4) 

Having computed as (^c), we can calculate the coupling constant at any scale in the 3- 
fiavor theory. For example, 

a(^) (/i = 2 GeV) = 0.2904. (A5) 



Appendix B: PERTURBATIVE RUNNING AND MATCHING FOR Zm 

In [l^ , the renormalization group equation for rrij-cn (/u) is solved to four- loop order 
(NNNLO) . Using the solution with our definition of the renormalization coefficients Z^, we 
obtain: 

rim /it) 

(Bl) 



^si ^ c {as (/ip) /tt) ^ri/mom ^^-^ 



c{as (/i) /vr) 



where the function c (x) is given by: 



c{x)={xy'{l+{-fi-Pao)x 
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with Pi = ^ and 
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The evaluation of the coefficients of the QCD (3 function and the running of are explained 
in Appendix |X] and the anomalous dimensions are 

RI/MOM(0) ^1 

Im 
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V 243 27 
where nj = 3. 

When applying Eq. (IBip . we need to choose a value of /iq, where the SI value is calculated. 
The exact value of fiQ is immaterial and for convenience we choose its value such that 



(a/io 



2. 



(B5) 



To match the renormalization coefficients _Zm from RI/MOM scheme to MS scheme, we 
have applied the three-loop matching factor [l4\ obtaining: 
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Appendix C: PERTURBATIVE RUNNING AND SCHEME MATCHING FOR Z„ 



The renormalization group equation for Zg is very similar to that for Zm [14] and we can 
reuse the solution of the equation from Appendix ([B]) to write: 
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where the function cl'^^l (x) has exactly the same functional form as c (x) defined in Eq. flB2p . 
but with the coefficients 7^ of the anomalous dimension 7^ replaced by those of 72: 
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where A^ = 3, which represents the number of colors, and nj = 3. 

When we match Zg from RI/MOM scheme to MS scheme, the three- loop matching factor 
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Appendix D: PERTURB ATI VE RUNNING AND SCHEME MATCHING FOR 

Zbk 



To remove (restore) the perturbative renorma 
one-loop renormalization group running formula 



ization group running of Zb^^, we use the 
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Appendix E: PERTURBATIVE RUNNING AND SCHEME MATCHING FOR Zt 



The anomalous dimension of the tensor current in the MS scheme is given at three-loops 
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For consistency we have adjusted the normahzation from that used in [27| such that 
satisfies the generic RG-equation for the renormahzation constant of the quark bihnear 

dlnZr , . 

= -7?^ a. - l?^al - 7?^a^ + 0{at) , (E2) 

with tts = tts/vr. 

The perturbative running for the tensor current has also been computed at three loops 



in the RI/MOM' scheme 271], and we use it to obtain the RI/MOM scheme anomalous 



dimension as follows. We consider the conversion function (7^^/^*-'^' ^ used to match the 



RI/MOM or RI/MOM' scheme to the MS scheme: 



^MS _ ^RI/MOMW^RI/MOMW ^j.g^ 

Applying the above renormahzation group equation (1E2I) we obtain 

o , ^ri/momC) 
^ri/momC) _ ^ms _ ^^^^r _ 

Since the only difi'erence between the RI/MOM and RI/MOM' schemes lies in the def- 
inition of the quark field renormahzation constants z^^^^*^^ and ^^^mom^ write 

^RI/MOMW ^ ^RI/MOMC) rp, , i. € • r , • • 

Cp = CrCg . ihe vertex part Cr 01 the conversion function is common 

to both the RI/MOM and RI/MOM' schemes. It follows that 

r.. ^RI/MOM 01 ^RI/MOM' 
RI/MOM' _ RI/MOM _ Cll'llL^2 111(^2 

'^^ (9 In /i^ d In /i^ 

RI/MOM RI/MOM' /-r^-x 

= 72 - 72 • (E5) 

Since both functions ^^^/^^^ an^j ^ri/mom known [14], we can now compute the anoma- 
lous dimension of the tensor current in the RI/MOM scheme from the known one in the 
RI/MOM' scheme. We note that since the r.h.s. of flESP is valid for any choice of F on the 
l.h.s., one may use the identity 

RI/MOM _ RI/MOM' / RI/MOM' RI/MOMx /-pp^x 

7r ~ 7r I7r' 7r' ) ■ \^^) 

T J . . RI/MOM 111 RI/MOM . ,17^ , RI/MOM' r 

in order to compute 7-^ here we have used 73 as m 11021) and 72 from 
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In this way we obtain the anomalous dimension: 
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from which we compute the running of Zt using ( 1B2[) . 

Combining (lA2p and (lE2p we compute the expression for the matching factor C^^' 
After expanding in we obtain: 
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Table I: The four factors Zs^p, Zy^j^, Zx and Zbi^ by which the matrix elements of the bare lattice 
bilinear operators and the ratio of matrix elements Bk should be multiplied in order to obtain 
the corresponding quantities renormalized in the RI/MOM or MS(NDR) schemes. The RI/MOM 
quantities are defined at a scale /x = 2.037 GeV, an available lattice momentum. The MS(NDR) 
quantities are provided at the scale /j, = 2 GeV. The first error given is statistical and the second 
systematic. This table summarizes the main results of this paper. 

Scheme Scale Zq Zs,p Zv,a Zt Zb^ 

RI/MOM 2.037 GeV 0.8086(28) (74) 0.466(14)(31) 0.7161(1) 0.8037(22) (55) 0.9121(38)(129) 
MS(NDR) 2.00 GeV 0.7726(30) (83) 0.604(18) (55) 0.7161(1) 0.7950(34) (150) 0.9276(52) (220) 



Table II: The factors, computed in perturbation theory, by which the matrix elements of the bare 
lattice operators should be multiplied in order to obtain those in the MS(NDR) scheme at the 
renormalization scale ji = 1.729 GeV. This table shows that the difference in the choice of the 
strong coupling constant leads to large uncertainty in the renormalization constants. 



Coupling 


Z5,p(1.729GeV) Zy,A(l-729GeV) Zr(l-729GeV) 


(1.729 GeV) 


aMF(1.729GeV) 


0.788 0.801 0.827 


0.979 


a^(1.729GeV) 


0.672 0.693 0.737 


0.963 



Table III: The perturbative renormalization constants at the conventional scale of = 2 GcV by 
renormalization group running from jj, = 1.729 GeV. The entries in the first column indicate which 
coupling was used in matching between the bare lattice operators and the MS(NDR) scheme at 
/i = 1.729 GeV. 

Coupling Z5,p(2GeV) Zv,A{2GeV) ZT(2GeV) ZB^(2GeV) 
aMF(1.729GeV) 0.822 0.801 0.813 0.993 

a^(1.729GeV) 0.701 0.693 0.725 0.977 
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Table IV: The quantity ]^Tr(S'[g^^^^) evaluated at the unitary mass points, rrivai = rrii and linearly 
extrapolated to the chiral limit mi = —rrires- 

{apf mi = 0.01 mi = 0.02 mi = 0.03 chiral limit 



0.347 0.0839(16) 0.1141(16) 0.1327(23) 0.0524(34) 

0.617 0.0558(13) 0.0810(15) 0.0980(20) 0.0283(28) 

0.810 0.0450(12) 0.0692(14) 0.0849(18) 0.0187(28) 

1.079 0.03744(82) 0.0583(11) 0.0741(16) 0.0130(20) 

1.234 0.0342(12) 0.0543(13) 0.0704(16) 0.0105(26) 

1.388 0.03203(84) 0.0512(11) 0.0665(15) 0.0092(20) 

1.542 0.03051(75) 0.04873(97) 0.0634(14) 0.0087(18) 

1.851 0.02640(92) 0.04472(99) 0.0597(13) 0.0047(20) 

2.005 0.02615(73) 0.04354(92) 0.0575(13) 0.0054(18) 

2.467 0.0236(10) 0.0104(10) 0.0510(13) 0.0010(20) 



Table V: The five bare vertex amplitudes Aj, z G {S,P,V,A,T} averaged over four sources, with 

m = mvai = 0.01. 

(ap)2 As Ap Ay Aa At 

0.347 2.125(86) 6.72(19) 1.1702(58) 1.0675(43) 0.8904(43) 
0.617 1.945(51) 4.45(11) 1.1419(37) 1.0938(30) 0.9404(26) 
0.810 1.856(37) 3.677(81) 1.1348(31) 1.1025(27) 0.9618(19) 
1.079 1.758(27) 3.022(57) 1.1335(29) 1.1135(27) 0.9882(16) 
1.234 1.715(24) 2.792(50) 1.1291(29) 1.1137(27) 0.9935(17) 
1.388 1.677(21) 2.600(43) 1.1328(26) 1.1191(24) 1.0065(13) 
1.542 1.642(19) 2.448(38) 1.1355(27) 1.1240(25) 1.0167(14) 
1.851 1.599(16) 2.239(32) 1.1387(29) 1.1301(27) 1.0310(16) 
2.005 1.578(15) 2.154(28) 1.1420(27) 1.1342(26) 1.0392(16) 
2.467 1.532(13) 1.979(23) 1.1495(29) 1.1434(29) 1.0577(19) 
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Table VI: The five bare vertex amplitudes Aj, i G {S*, P, V,^, T} averaged over four sources, with 
mi = ruyax = 0.02. 



As 



Ap 



At 



0.347 1 

0.617 1 

0.810 1 

1.079 1 

1.234 1 

1.388 1 

1.542 1 

1.851 1 

2.005 1 

2.467 1. 



.828(45 
.774(30 
.721(24 
.655(19 
.637(16 
.608(15 
.581(14 
.552(11 
.532(11 
1981(91) 



5.09(14) 

3.600(82) 
3.052(61) 
2.590(45) 
2.428(40) 
2.283(33) 
2.175(30) 
2.019(24) 
1.955(23) 
1.819(18) 



1.1745(46) 

1.1465(30) 
1.1360(24) 
1.1331(22) 
1.1307(21) 
1.1323(21) 
1.1351(21) 
1.1389(22) 
1.1416(23) 
1.1198(26) 



1.0412(28) 
1.0838(21) 
1.0943(19) 
1.1069(20) 
1.1083(20) 
1.1141(19) 
1.1199(20) 
1.1275(21) 
1.1315(22) 
1.1122(25) 



0.8930(31) 

0.9414(18) 
0.9614(15) 
0.9870(12) 
0.9930(12) 
1.0049(11) 
1.0159(12) 
1.0309(12) 
1.0390(14) 
1.0580(17) 



Table VII: The five bare vertex amplitudes Aj, i G {S, P, V, A, T} averaged over four sources, with 
mi = mvai = 0.03. 

(ap)2 As Ap Av A A At 



0.347 1.723(56) 4.10(14) 1.1809(56) 1.0357(23) 0.9020(24) 

0.617 1.702(37) 3.049(87) 1.1457(37) 1.0769(19) 0.9451(17) 

0.810 1.663(28) 2.658(66) 1.1356(31) 1.0886(17) 0.9642(14) 

1.079 1.610(21) 2.307(49) 1.1325(27) 1.1015(18) 0.9883(12) 

1.234 1.591(18) 2.182(42) 1.1294(25) 1.1050(20) 0.9951(12) 

1.388 1.569(15) 2.076(37) 1.1312(25) 1.1105(20) 1.0061(11) 

1.542 1.548(13) 1.991(33) 1.1337(26) 1.1157(21) 1.0161(12) 

1.851 1.520(10) 1.869(27) 1.1366(26) 1.1228(22) 1.0300(13) 

2.005 1.5065(96) 1.820(25) 1.1395(27) 1.1271(24) 1.0382(15) 

2.467 1.4764(78) 1.717(20) 1.1464(27) 1.1371(26) 1.0561(17) 
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Table VIII: Groups of non-exceptional momenta satisfying Pi = P2 = {Pi ~ P2)^- The individual 
integers {rix, riy, Uz, nt) should be multiplied by 2n/L(i, with = Ly = Lz = 16 and Lt = 32. 
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l^-i,-i,-Z,U j 


(-l,U,-l,-4j 




( ^ 1 9 ^^ 

1-1,-1, -z,u 1 


('in 1 /i^ 




( ■[ 1 9 n^ 


("n 1 1 /l^ 




(.^ _i -9 

V 5 5 5 / 


(T) -1 -1 4"! 




(-1,-1,-2,0) 


(0,1,-2,-2) 




(-1,-1,-2,0) 


(0,1,-2,2) 




(-1,-1,-2,0) 


(1,-2,-1,0) 




(-1,-1,-2,0) 


(1,0,-2,-2) 




(-1,-1,-2,0) 


(1,0,-2,2) 
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Table IX: Groups of non-exceptional momenta satisfying Pi = P2 = {pi ~P2)'^, continuing 
ble lvmi 



(apf 


P\ 


P2 


1.234 


(0,2,2,0) 


(2,2,0,0) 




(0,2,2,0) 


(0,2,0,4) 

V 7 7 7 / 




(0 2 2 0) 


(0 2 4) 




(0,2,2,0) 


(-2,2,0,0) 




(0,2,2,0) 


(0,2,0,-4) 




(0,2,2,0) 


(2,0,2,0) 




(0,2,2,0) 


(0,0,2,-4) 




(0,2,2,0) 


(-2,0,2,0) 


1.542 


(1,1,2,4) 


(2,1,2,-2) 




(1,1,2,4) 


(1,-2,2,2) 




(112 4) 


f-2 1 2 21 




(1,1,2,4) 


(-2,1,1,4) 




(1,1,2,4) 


(1,2,2,-2) 




(1,1,2,4) 

\ 7 7 7 / 


(1,-2,1,4) 




(1,1,2,4) 


(2,1,-1,4) 




(1,1,2,4) 


(1,2,-1,4) 


2.467 


(2,2,2,4) 


(2,2,-2,4) 




(2,2,2,4) 


(2,-2,2,4) 




(2,2,2,4) 


(-2,2,2,4) 




(2,2,2,4) 


(2,2,2,-4) 
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Table X: Results from fitting the coefficient for mass term in (A^ — Ay)/[(A/i + Ay)/2]. The linear 
dependence is assumed to be ci and the quadratic dependence is assumed to be C2^. The 

respective ^ /d.o.f is also listed. Both the coefficient ci more nearly agreeing with its expected 
value of 1 and the smaller suggest that the linear description is to be preferred. We use the 
value Aqcd = 319.5 MeV. 







VA / / i 


Co 




0.347 


-3.84(75) 


3.0(3.4) 


-14.7(3.1) 


6.6(5.1) 


0.617 


-3.33(67) 


2.2(2.8) 


-12.9(2.8) 


5.4(4.4) 


0.810 


-3.06(56) 


1.2(2.1) 


-12.1(2.5) 


4.1(3.7) 


1.079 


-3.01(42) 


0.4(1.3) 


-12.4(2.0) 


3.3(3.6) 


1.234 


-2.96(47) 


6.2(5.0) 


-11.4(2.1) 


12.8(7.2) 


1.388 


-2.58(36) 


1.7(2.7) 


-10.5(1.7) 


6.2(4.8) 


1.542 


-2.49(34) 


0.4(1.4) 


-10.2(1.6) 


3.3(3.6) 


1.851 


-2.33(35) 


0.06(41) 


-9.5(1.6) 


1.8(2.4) 


2.005 


-2.21(28) 


0.02(23) 


-9.2(1.3) 


1.3(2.2) 


2.467 


-1.89(32) 


0.01(22) 


-7.4(1.4) 


0.7(1.5) 
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Table XI: Values for |(Ayi + Ay), A5, and At extrapolated to the chiral limit using a linear mass 
fit. 



(ap)2 i(AA + Av) 



A5 



At 



0.347 1.1211(56 

0.617 1.1226(49 

0.810 1.1228(41 

1.079 1.1275(43 

1.234 1.1242(44 

1.388 1.1292(40 

1.542 1.1333(43 

1.851 1.1381(47 

2.005 1.1417(46 

2.467 1.1504(51 



2.28(14) 0.8800(66) 
2.060(88) 0.9363(36) 
1.952(66) 0.9593(25) 
1.836(47) 0.9873(22) 
1.784(43) 0.9915(22) 
1.736(37) 1.0061(18) 
1.694(32) 1.0169(20) 
1.644(27) 1.0319(24) 
1.617(25) 1.0400(27) 
1.564(22) 1.0593(33) 
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Table XII: The non-perturbative factor Z^^^^^*^^^ as a function of the scale calculated from 
A5 and the corresponding values for Z^^. Note that the values for given in column three 
are obtained from those in column two by applying the RI/MOM — MS perturbative matching 
factors after the O^ajj)'^ lattice artifacts have been removed using an intermediate conversion to a 

scale-invariant scheme as described in the text. 



/i(GeV) Z^^/^O^ Z^S 

1.018 2.85(18) 1.625(47) 

1.358 2.56(11) 1.758(51) 

1.556 2.428(80) 1.731(51) 

1.796 2.273(56) 1.690(49) 

1.920 2.216(49) 1.669(49) 

2.037 2.146(42) 1.651(48) 

2.147 2.087(37) 1.634(48) 

2.352 2.018(29) 1.605(47) 

2.448 1.978(27) 1.593(46) 

2.716 1.899(22) 1.562(46) 
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Table XIII: The non-perturbative factor Zq as a function of the scale /i calculated from 

^ (A^ + Ay) and the corresponding values for Z^^. Note that the values for Z^^ given in column 
three are obtained from those in column two by applying the RI/MOM — MS perturbative matching 
factors after the O^ajj)'^ lattice artifacts have been removed using an intermediate conversion to a 
scale-invariant scheme as described in the text. 

^(GcV) zf/^o^ Zf^ 

1.018 0.8028(40) 0.8010(31) 

1.358 0.8039(35) 0.7849(30) 

1.556 0.8041(30) 0.7798(30) 

1.796 0.8074(31) 0.7754(30) 

1.920 0.8050(32) 0.7736(30) 

2.037 0.8086(28) 0.7722(30) 

2.147 0.8115(31) 0.7710(30) 

2.352 0.8150(34) 0.7691(29) 

2.448 0.8176(33) 0.7684(29) 

2.716 0.8238(37) 0.7665(29) 
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Table XIV: The non-perturbative factor function of the scale \i calculated from 

Ay and the corresponding values for Note that the values for given in column three 

are obtained from those in column two by applying the RI/MOM — MS perturbative matching 
factors after the O^ajj)'^ lattice artifacts have been removed using an intermediate conversion to a 
scale-invariant scheme as described in the text. 

^(GcV) zf'^'''^ 

1.018 0.9121(74) 0.8812(38) 

1.358 0.8583(46) 0.8355(36) 

1.556 0.8380(32) 0.8194(35) 

1.796 0.8177(27) 0.8048(34) 

1.920 0.8118(27) 0.7986(34) 

2.037 0.8037(22) 0.7935(34) 

2.147 0.7981(21) 0.7892(34) 

2.352 0.7899(18) 0.7821(33) 

2.448 0.7862(17) 0.7791(33) 

2.716 0.7779(16) 0.7719(33) 



Table XV: The quantity evaluated at the unitary points where myai = mi = m. 

fi{GeY) m = 0.01 m = 0.02 m = 0.03 

0.954 0.9663(69) 0.9737(52) 0.9538(44) 

1.272 0.9347(39) 0.9387(35) 0.9315(30) 

1.458 0.9266(31) 0.9289(30) 0.9245(26) 

1.683 0.9189(25) 0.9189(25) 0.9167(24) 

1.799 0.9151(22) 0.9137(23) 0.9126(23) 

1.909 0.9114(23) 0.9106(22) 0.9102(22) 

2.012 0.9085(22) 0.9077(20) 0.9078(21) 

2.204 0.9045(23) 0.9026(20) 0.9035(21) 

2.294 0.9018(20) 0.9004(19) 0.9020(20) 

2.545 0.8974(21) 0.8953(19) 0.8978(19) 
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Table XVI: The non-perturbative factor ^^j^^^M ^ function of the scale /j, and the corresponding 
values for Z^^. Note that the values for given in column three are obtained from those 
in column two by applying the RI/MOM — MS perturbative matching factors after the 0{an)^ 
lattice artifacts have been removed using an intermediate conversion to a scale-invariant scheme 

as described in the text. 

MGeV) Zg^/^Q^ 

1.018 0.985(11) 1.0016(56) 

1.358 0.9397(61) 0.9651(54) 

1.556 0.9295(48) 0.9507(54) 

1.796 0.9208(42) 0.9370(53) 

1.920 0.9168(38) 0.9311(52) 

2.037 0.9121(38) 0.9261(52) 

2.147 0.9088(37) 0.9217(52) 

2.352 0.9045(39) 0.9145(52) 

2.448 0.9011(35) 0.9114(51) 

2.716 0.8961(38) 0.9038(51) 
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Figure 1: The quantity j^Tr (5'jj^t\) plotted versus (ap)^ for the unitary mass points mi = 0.01, 
0.02 and 0.03 and at the hnearly extrapolated, chiral Hmit mi 
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Figure 2: The ratio (jVA+Ay^ /2 P^o^^^*^ ^ ^ function of momentum at the unitary mass points 
n^-vai = TTii and in the chiral hmit evaluated by linear extrapolation in mi. The 5-10% difference 
at low momentum decreases rapidly as the momentum increases. At the scale ~ 2 GeV, or 



(ap)^ ~ 1.4, the difference is about 1%, which contributes to the systematic error in Zbj. 
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Figure 3: The difference A^i — Ay computed using four different quenched DBW2 lattice ensembles. 
These ensembles have quite different lattice scales. In addition the values of Lg, the extent in the 
5th dimension used in computing the DWF propagators, also varies significantly. This provides 
compelling evidence that the observed chiral symmetry breaking is not an explicit breaking from 
finite Ls, but rather represents the high energy tail of QCD dynamical chiral symmetry breaking 
which would vanish if we were able to p erform the NPR calculation at high enough energy. The 



data shown come from Refs. 



151, 



281, 



29| 
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Figure 4: The x^/d.o.f which results from fitting the momentum dependence of the quantity 
A ^-Ay — (extrapolated to the chiral limit) to the form p^". We conclude that the best choice for 



(A.4+Av)/2 

n lies between 2 and 3 and that it is unlikely that the term {qq)"^ /p^ gives the dominant contribution 
to this chiral symmetry breaking. 
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Figure 5: The division of a general vertex graph into subgraphs. If the four-legged, internal 
subgraph r2 carries momenta p ~ Aqcd it can introduce low energy, (8, 8) chiral symmetry breaking 
into such an amplitude even in the limit that the momenta external to the entire diagram F, 

included in the outer dashed box, grow large. As discussed in the text, such a limit will be 
suppressed by if the external momenta are non-exceptional but by only l/p^ for the exceptional 
case. 




Figure 6: Sample diagram in which two low-momentum {k 2± Aqcd) fermion propagators appear 
in a graph which is suppressed at high momentum only by a single factor oil/p^. 
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Figure 7: The value of 



A4-A1, 



calculated with non-exceptional kinematics, which requires the 



(Aa+Av)/2 

sum of any subset of external momenta be non-zero. With this condition the chiral symmetry 
breaking is highly suppressed (as compared to Fig. [2]) and vanishes almost completely over the 
available momentum region. 
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Figure 8: The average ^ {Aa + Ay) plotted as a function of momentum and evaluated for a unitary 
choice of masses and in the chiral limit. The chiral limit is taken using a linear fit. 
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Figure 9: Comparison of linear (eq. (f39|) - top panel) and quadratic (eq. (f40]) - bottom panel) fits to 
the dependence of the chiral symmetry breaking difference (Aa+AvO/2 °^ quark mass ni^oi = mi 
at the scale /U = 2.04 GeV. These plots suggest that a linear description is more accurate. This 
conclusion is borne out by the properties of the actual fits shown in Table |Xl 
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Figure 10: A plot showing the linear extrapolation of ^ {Aa + Ay) (evaluated at the scale fj, 
2.04 GeV) to the chiral limit. The three data points are evaluated at the unitary points m-vai = " 
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Figure 11: The difference between the quantities Ap and -n, ^ V , divided by their average 

■•■ iZ TTi; -|-7Ti]-cs 

is plotted versus momentum for unitary quark masses. This provides a test of the axial Ward- 
Takahashi identity. 
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Figure 12: The difference between the quantities A5 and j2~~^n^Y^^ divided by their average 
for each sea quark mass. The difference appears to zero within errors. This is a test of the vector 



Ward-Takahashi identity. The plot uses propagators from three sources. 
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Figure 14: Momentum dependence of the double pole coefficient, Cdp, fit to the expected 
behavior. Good agreement is seen. 
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Figure 15: The quantities Z^^^'^^^ (/i) and (fi) plotted versus the square of the scale afi. 
Here Z^ (/i) is obtained by dividing Z^^^*^^ (/x) by the predicted perturbative running factor. 
Shown also is the linear extrapolation of Z^ (fj,) = Z^ + c (fl/u)^ using the momentum region 
1.3 < {a^ly < 2.5 to remove lattice artifacts. 
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Figure 16: The mass renormalization factor expressed in the MS scheme. These results are 
obtained by applying the perturbative running factor to Z"^. The value we are interested in is 
Z^ (/X = 2 GeV). The upper and lower curves show the statistical errors. 
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Figure 17: The quantities ^^V^^OM ^gj plotted versus the square of the scale afi. 

Here Z^^ {^) is obtained by dividing ^^^mom predicted perturbative running factor. 

Shown also is the linear extrapolation of Z^^ (fj,) = Z^^ + c (a^u)^ using the momentum region 
1.3 < {a^ly < 2.5 to remove lattice artifacts. 
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Figure 18: The wave function renormalization factor Zq expressed in the MS scheme. These results 
are obtained by applying the perturbative running factor to Z^^. The value we are interested in is 
Z^'^ (/X = 2 GeV). The upper and lower curves show the statistical errors. 



66 




67 



0.9 



N 



0.8 



0.7 



C) 



o 


RI/MOM 


□ 


SI 




(a\i f fit 



□gn rff an nm m ^ m 



0.5 



1.5 



(afx) 



2.5 



Figure 20: The quantities ^^^Z^*-*^ and (fi) plotted versus the square of the scale a/i. 
Here Zj^ {fi) is obtained by dividing ^^^/'^'-'^ (^^^ ]-,y i^j^e predicted perturbative running factor. 
Shown also is the linear extrapolation of Zj^ (/x) = Zj^ + c (afi)'^ using the momentum region 



1.3 < {a^ly < 2.5 to remove lattice artifacts. 
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Figure 21: The wave function renormalization factor expressed in the MS scheme. These results 
are obtained by applying the perturbative running factor to Z"^ . The value we are interested in is 
Z^"^ = 2 GeV). The upper and lower curves show the statistical errors. 
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Figure 22: The mixing coefficient FMyy^^^ vv-AA ^'^'^ ^^'^ three unitary mass values and linearly 
extrapolated to the chiral limit. 
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Figure 24: A possible identification of subgraphs appearing in the chirahty violating mixing between 
Oil and other four-quark operators. The disconnected subdiagram V has degree of divergence 
d = —2 for the case of exceptional momenta shown here. This permits a complex pattern of low- 
energy, vacuum chiral symmetry breaking coming from the low-energy, four-quark subgraph Fi to 
enter such an amplitude with only a mild suppression. 
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Figure 25: The mixing coefficient FMyy^j^^ vv-AA calculated at non-exceptional momenta. 
When extrapolated to the chiral limit the mixing coefficient vanishes, which shows that chiral 
symmetry breaking as shown in Fig. [22] comes from the existence of a low-energy sub-diagram that 
enters because of the special choice of external momenta. 
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Figure 26: The mixing coefficient FMyy_^_^^ ss-pp unitary choices of the mass. 
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Figure 27: The mixing coefficient FMyy_^j^j^ ss-pp calculated at non-exceptional momenta. When 
extrapolated to the chiral limit the mixing coefficient vanishes, which shows that chiral symmetry 
breaking as shown in Fig. [26] comes from the existence of a low-energy sub-diagram that enters 
because of the special choice of external momenta. 
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Figure 28: The mixing coefficient FMyy^^^ ss+pp ^'-'^ unitary choices of the mass. The coefficients 
are very tiny over the region of medium to large momenta. 
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Figure 29: The mixing coefficient FMyy^^^ rpj, for unitary choices of the mass. These coefficients 
agree weh with zero. 
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Figure 31: The renormalization factor ^^^/^'-'^ evaluated for unitary mass values and extrapolated 



to the chiral limit. 
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Figure 32: The quantities _^^Vi^OM plotted versus the square of the scale a/x. Here Zf^ 

is obtained by dividing z^^^^*^^ ]-,y ^j^g predicted perturbative running factor. Shown also is the 
linear extrapolation of = + c{a^) using the momentum region 1.3 < (anY < 2.5 to 

remove lattice artifacts. 
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Figure 33: The renormalization factor Zbj^ expressed in the MS scheme. These results are ob- 
tained by applying the perturbative running factor to Z^^. The value we are interested in is 
= 2 GeV). The upper and lower curves show the statistical errors. 
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